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Abstract 

The uncertainty in the absolute value of the fluorescence yield is 
still one of the main contributions to the total error in the reconstruc- 
tion of the primary energy of ultra-energetic air showers using the 
fluorescence technique. A significant number of experimental values 
of the fluorescence yield have been published in the last years, however 
reported results are given very often in different units (photons/MeV 
or photons/m) and for different wavelength intervals. In this work 
we present a comparison of available results normalized to its value 
in photons/MeV for the 337 nm band at 800 hPa and 293 K. The 
conversion of photons/m to photons/MeV requires an accurate deter- 
mination of the energy deposited by the electrons in the field of view 
of the experimental setup. We have calculated the energy deposition 
for each experiment by means of a detailed Monte Carlo simulation 
including when possible the geometrical details of the particular setup. 
Our predictions on deposited energy, as well as on some geometrical 
factors, have been compared with those reported by the authors of the 
corresponding experiments and possible corrections to the fluorescence 
yields are proposed. 



1 Introduction 



Charged particles of a cosmic ray shower, mostly electrons, passing through 
the atmosphere lose energy by inelastic collisions with air molecules. A small 
fraction of the total deposited energy is emitted by molecular nitrogen as 
UV fluorescence radiation in the ~ 300 — 400 nm spectral range. The de- 
tection of this radiation provides a precise determination of the longitudinal 
profile allowing the reconstruction of the cosmic ray properties, in particular 
its energy. This technique [T] was first successfully used by the Fly's Eye 
experiment [2] and later by HiRes [3J. Fluorescence telescopes are presently 
being employed by the Pierre Auger Observatory [3] and the Telescope Ar- 
ray experiment [5]. On the other hand, the JEM-EUSO project [6] and the 
S-EUSO free-flying satellite mission [7j are being designed for the detection 
of fluorescence traces of air showers viewing downward from the top of the 
atmosphere. 

The ratio between number of fluorescence photons and deposited energy, 
i.e., the fluorescence yield Y, is a key calibration parameter determining the 
energy scale of the fluorescence telescopes. A number of absolute measure- 
ments of the fluorescence yield have been carried out in laboratory experi- 
ments in the last years. 

At atmospheric pressure, fluorescence emission in the spectral range of 
interest basically comes from the Second Positive (2P) system of N2 and, 
to a much less extent, the First Negative (IN) system of Nj ■ Each one 
results from the transition between two electronic states and consists of a set 
of molecular bands corresponding to different combinations of vibrational 
levels v-v' of the respective upper and lower states. Excitation cross section 
of the IN system decreases with energy in the 10 2 — 10 6 eV range showing 
a very smooth growing behavior at larger energies. That of the 2P system 
peaks at about 15 eV decreasing strongly with an E~ 2 dependence. As 
a consequence, high-energy electrons themselves are very inefficient for the 
generation of air fluorescence. In fact, as is well known, fluorescence emission 
along the track of an energetic primary electron is mainly induced by low- 
energy secondary electrons ejected in successive ionization processes. The 
efficiency for fluorescence emission varies with the energy spectrum of these 
secondarie^J which in principle depends on the primary energy. However, a 

1 For instance, electrons with energies below the threshold for fluorescence excitation 
promptly lose all their energy, but induce no fluorescence. 
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theoretical calculation using a detailed Monte Carlo simulation [HI El [TO] has 
demonstrated that the fluorescence yield is nearly constant for energies above 
a few kiloelectronvolts. This E- independent behavior of the fluorescence yield 
has been checked experimentally [11] [131 EG] • 

A number of experimental values of the absolute fluorescence yield are 
available in the literature. Essentially, in these experiments a beam of elec- 
trons crosses a collision chamber filled with air at known conditions. However 
different techniques can be used. Many authors [131 [HI [151 EE El HH] have 
used electrons from a source of 90 Sr with energy around 1 MeV. Other abso- 
lute measurements have been performed with higher energy electrons from 
accelerators. The MACFLY Collaboration [T3] used the CERN/SPS-X5 test 
beam facility which delivers a pulsed electron beam of about 10 4 electrons per 
spill (4.8 s duration) every 16.8 s. Measurements at 20 and 50 GeV were car- 
ried out using this facility. The FLASH Collaboration [191 120] used the Final 
Focus Test Beam facility at the Stanford Linear Accelerator Center which 
provided 28.5 GeV electrons in 3 ps pulses of about 10 8 electrons at a rate of 
10 Hz. The AIRFLY Collaboration is carrying out absolute measurements 
of the air-fluorescence yield using different accelerators. Preliminary results 
obtained with 350 MeV electrons at the Beam Test Facility of the Instituto 
Nazionale di Fisica Nucleare have been presented in [21]. New measurements 
using a 120 GeV proton beam from the Meson Beam Test Facility at Fermi- 
lab were presented at the 6th AFW [22J and very likely final results will be 
published soon. 

In this work we will compare the absolute values reported by the above 
authors. In many occasions comparison cannot be done directly. For in- 
stance, some authors measure single intense fluorescence bands while others 
detect the integrated fluorescence in a wide spectral range. Comparison be- 
tween these experimental results can be carried out if the relative intensities 
of the bands along the fluorescence spectrum are known. We will normalize 
all the results to the most intense fluorescence band at 337 nm using ac- 
curate experimental data in full agreement with theoretical intensity ratios 
proposed in El ED] • 

On the other hand, some authors report their measurements in units of 
photons per electron and meter while others give directly the Y parameter 
in photons per unit deposited energy (e.g., photons/Me V). The conversion 
factor between both magnitudes is the energy deposited per electron and 
unit path length. As it has been pointed out [SI |H1 EE] , an accurate deter- 
mination of the deposited energy requires discounting that corresponding to 
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secondary electrons outside the field of view of the optical system used to 
detect the emitted fluorescence. To evaluate this effect we have carried out 
a MC simulation of the various experimental setups. When possible a dedi- 
cated simulation including geometrical details has been performed. Also the 
effect of the spatial distribution of deposited energy on the optical efficiency 
of the experimental setup has been evaluated in some cases. As a result of 
this study absolute values of the fluorescence yield normalized to that of the 
337 nm band at 293 K and 800 hPa in dry air will be shown with a discussion 
on possible corrections to be applied. A previous comparison at 1013 hPa 
for some of the above measurements has been published recently [10J. In 
this work we will update those results and extend the comparison to other 
experimental results. 

2 Fluorescence Yield 
2.1 Units 

In the literature, several parameters have been used to measure the amount 
of fluorescence light induced by an electron moving in the atmosphere. A 
detailed description of the various magnitudes can be found in [23]. The 
number of photons emitted per electron and unit path length have been 
used extensively. We will name this parameter as e and, for a band v — v', 
it is proportional to the number of nitrogen molecules per unit volume N 
and the reciprocal of the 1 + P/ P' v Stern- Volmer factor, accounting for the 
collisional de-excitation of the upper level v, where P' v is a characteristic 
pressure depending on the gas composition and temperature. If ignoring the 
contribution of secondary electrons, the proportionality constant would be 
the so-called optical cross section. 

While most secondary electrons have low energy, and thus, short range at 
atmospheric pressure, a small number of high-energy secondaries can escape 
the observation region. As a consequence, in a laboratory experiment a frac- 
tion of the fluorescence light cannot be detected. Following the parametriza- 
tion in [8] [HI EI] , an effective optical cross section cr°^, can be defined in such 
a way that the total number of photons of the v — v' band detected under 
given experimental conditions (including contributions from both primary 
and secondary electrons) can be expressed as 
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of 



£vv , = N ^— . (1) 

1 + PjP' v 1 ; 

The parameter a^, depends on E, P and the geometrical features of 
the experimental setup. Therefore, even for the same primary energy, the 
measured e m > value depends on the particular experiment. In addition, the 
pressure dependence of would be predicted by the Stern- Volmer law 

is distorted in a geometry-dependent way. 

A further pressure dependence arises from the vibrational relaxation of 
excited nitrogen induced by collisions with surrounding nitrogen molecules in 
the ground state [21] . The effect of vibrational relaxation is twofold. Firstly, 
it contributes to the total quenching of a given level v and thus is included in 
the Stern- Volmer factor ([I]). Secondly, vibrational relaxation of upper levels 
(> v) provides an additional excitation channel. Although, in principle, 
this effect may lead to a departure from the Stern- Volmer law, significant 
deviations have not been appreciated experimentally for any N2 band in 
air [25]. In fact, for the 2P(0,0) band, for which detailed data on vibrational 
quenching are available [24j . a simple calculation shows that the only relevant 
effect of vibrational quenching is a variation of the resulting characteristic 
pressure respect to the one from pure electronic quenching |2"5] . 

The number of photons per unit energy deposited in the medium is a more 
appropriate parameter for use in fluorescence detectors of cosmic ray showers. 
Following an already common usage, we will keep the term fluorescence yield 
to refer to this magnitude Y. For a given band, e vv i and Y vv i are related by 

Y , - — (2) 

where (dE/dx)d ep is the energy deposited per primary electron and unit path 
length. In the above equation, both the emitted photons and the deposited 
energy should correspond to the same volume, which is defined by the field of 
view of the experimental detection system. Note that a significant fraction of 
energy is deposited by secondary electrons which move away from the path 
of the primary electron. It has been shown [SI HO] that Na^, and (dE/dx)d ep 
are proportional, resulting Y vv > to be given by 

Y°, 

Yvv ' = 1 + pik* ' (3) 
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where the Y® v , parameter is independent of the experimental conditions and 
represents the fluorescence yield at null pressure, i.e., in the absence of 
quenching. 

While measurements of the fluorescence yield at atmospheric pressure by 
different authors show a reasonable agreement, the available data on the 
parameters show large discrepancies Therefore the extrapolation of the 
fluorescence yield to null pressure is unsafe, making the Y® v , parameter un- 
suitable for our proposes. In this work, comparison will be thus performed 
using the fluorescence yields measured at given atmospheric conditions. 

2.2 Deposited energy 

The energy deposited by an electron per unit path length within a given 
volume can be expressed [10] as 



(S ) d = n * ^ d ° ep) + (jBdcp) ] ai ° n ' (4) 

where iV air is the density of air molecules and <j- lon is the ionization cross 
section, including the density correction for high energies. The parameters 
(-Eyep) and (E^ep) are the average values of the energy deposited per ioniza- 
tion process within the observation volume by the primary electron itself and 
the secondaries, respectively. For electron energies above 410 eV, contribu- 
tion from X-rays arising from K-shell ionization should be also added to the 
energy deposition [lOj. 

The (E® ) parameter can be calculated from 

(^ ep) = (^ X c)^ + ra+/, (5) 

^ion 

where (E™) and / are the average excitation energy and the ionization po- 
tential of the ionized molecule, respectively, both contributing to the energy 
deposition of the primary electron after an ionizing interaction. The first term 
in the right hand of the equation accounts for the contribution from direct 
excitations without ionization, being (E exc ) the average excitation energy of 
the molecule and <r exc the excitation cross section for a primary electron. 



2 Notice that many authors ignore the additional pressure dependence of e vv > due to 
secondary electrons [TP] . 
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The parameters involved in the above expressions are either available in 
the literature or can be inferred from molecular data, except for (E^cp), which 
has to be calculated from a simulation of the secondary processes inside the 
observation volume. 

Next we will describe a Monte Carlo algorithm to calculate both the en- 
ergy deposited by secondary electrons and their contribution to the fluores- 
cence emission for a simple geometry. In section [3] a more detailed simulation 
including the geometrical features of several experiments will be presented. 

2.3 Monte Carlo analysis under simple geometrical 
assumptions 

Ignoring fine geometrical details, the observation region of an experimen- 
tal setup can be described by a sphere of radius R. In addition, we will 
assume that primary electrons interact at its center from where secondaries 
are ejected and tracked until they escape the sphere or are stopped. A Monte 
Carlo algorithm for the calculation of both deposited energy and fluorescence 
emission in this simple geometry has been described in detail in [HI OH] • Next 
the main features and ingredients of the simulation will be outlined and some 
improvements described. 

The algorithm generates secondary electrons with energies distributed as 
given by the extended Opal formula proposed in [10J, which depends on the 
input primary energy. Each generated electron takes a step of random size 
with mean value {N^a) -1 , where a is the total cross section for the electron 
energy. After this step, the type of process is randomly chosen between elas- 
tic, excitation, ionization and bremsstrahlung, being their relative weights 
the corresponding cross sections. While there are no energy losses in elastic 
interactions, in excitation and ionization processes the electron is assumed to 
locally deposit (E exc ) and (E™) + /, respectively. In case of ionization, the 
electron loses an additional amount of energy to generate a new secondary, 
whose energy is randomly calculated from the above-mentioned extended 
Opal formula. With a probability corresponding to the cross section for K- 
shell ionization, also a 410 eV X-ray can be generated at the interaction point, 
producing another ionization at a random distance^] of mean value given by 

3 In previous versions of our Monte Carlo algorithm [10] , X-rays were not propagated, 
assuming all their energy to be locally deposited without producing further ionizations. 
While this approximation has no significant effect on the deposited energy, due to the short 
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the photon attenuation coefficient from [27J . In case of bremsstrahlung, the 
energy loss is estimated as 3% of the electron energy and the produced pho- 
tons are assumed to escape the sphere^] If an excitation (ionization) process 
occurs, the amount of 2P (IN) fluorescence light emitted is calculated from 
the ratio of the 2P (IN) cross section to the total excitation (ionization) 
one. The direction of the electron after an interaction is calculated from an 
approximated differential cross section [2H]. The stepped trajectory of each 
electron is so iteratively simulated until it escapes the sphere or its energy 
is below 11 eV (i.e., the threshold for fluorescence production), when the 
electron is considered to be stopped and deposits its remainder energy. New 
secondary electrons arising from ionization processes are added to the list of 
electrons to be simulated following the same procedure. 

Input parameters in the algorithm are the energy of the primary electron, 
the density of air molecules (or pressure at given temperature) and the radius 
of the sphere. As a result of the simulation the value of (E dep ) is obtained, 
and then, the total energy deposited per primary path length is calculated 
from expressions (Jij) and (J5|. Likewise, the a^, parameter can be calculated 
from the average number of fluorescence photons emitted inside the sphere 
in the absence of quenching resulting from the simulation. As expected, 
both magnitudes turn out to be a function of the product R ■ N aiT , that is, 
proportional to the number of molecules that a secondary electron encounters 
in its way out of the sphere, and the quotient provides the Y® v , parameter, 
which is basically independent of E and R ■ N^. 

As discussed in [10], our calculation of the fluorescence emission is sub- 
jected to a systematic uncertainty, that we estimate in about a 25%, due to 
the uncertainties in the various molecular parameters. However the evalua- 
tion of the energy deposited in the sphere is much more accurate. In fact, 
in the R ■ A^ air — > oo limit, for which no secondary escapes, the predictions 
from the simulation are in perfect agreement with the Bethe-Bloch formula 
for the collisional energy loss [29J including the corresponding density cor- 
rection. For this result, a crucial ingredient in the simulation is the energy 
distribution of secondary electrons. The extended Opal formula proposed 
in [10] provides this spectrum for a very wide energy interval (eV - GeV) of 

range of the X-rays, it leads to an underestimation of about 10% for the total fluorescence 
emission. 

4 We checked that contribution from bremsstrahlung photons to the energy deposition 
is negligible even for large experimental setups and GeV electrons, where energy losses 
due to bremsstrahlung are dominant. 
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both primary and secondary electrons. 

Simulation results of the deposited energy for several R ■ values have 
been presented in [10]. For the purpose of comparison, in figure [I] the energy 
deposited by an electron per unit mass thickness^ in dry air at atmospheric 
pressure is represented against the primary energy for R = 1 cm and 10 cm. 
Below 0.5 MeV the fraction of energy deposited outside the sphere is negli- 
gible at this pressure and therefore the energy deposition for both R values 
approaches the energy loss. As can be seen in the figure, the deposited en- 
ergy is very weakly dependent on the size of the observation region at given 
pressure. For instance, the energy deposited by electrons of 1 MeV (1 GeV) 
increases by 9% (6%) when R increases from 1 to 10 cm at atmospheric pres- 
sure. This justifies that a reasonable accuracy in the deposited energy can be 
achieved even though the fine geometrical details of the experimental setup 
are not taken into account. 

Notice that in our algorithm, individual electron-molecule collisions as 
well as the emission of IN and 2P fluorescence are simulated with a cutoff 
energy of 11 eV. On the contrary, usual MC codes (e.g., GEANT and EGS) 
transports particles in macroscopic steps adding up the energy deposition 
in each step. Fluorescence emission, which is mostly generated by electrons 
with energies of about a few tens of electronvolts, cannot be calculated by 
these codes unless the user assumes a certain value of the fluorescence yield. 

2.3.1 Comparison of results on deposited energy with other 
simulations 

The experimental determination of the fluorescence yield requires the calcu- 
lation of the energy deposited in the observation volume. Consequently, any 
possible error in the evaluation of this magnitude translates to the fluores- 
cence yield. Several experiments have provided results on their simulations 
which can be compared with our predictions for typical sizes of experimental 
setups (see figure [T]). 

In their paper of 2006 the FLASH Collaboration [19] showed results on 
deposited energy versus E for a 1 cm thick slab of air at atmospheric pressure 
obtained with GEANT3, which are significantly larger than our predictions 
for R = 1 cm. In a more recent paper [20], they obtained an energy deposition 

Results expressed in MeV per unit mass thickness are valid for both dry air and pure 
nitrogen since inelastic cross sections are approximately proportional to the atomic number 
Z, which is a factor of two lower than the mass number A for both Nitrogen and Oxygen. 
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Figure 1: Energy deposited per electron and unit mass thickness at atmo- 
spheric pressure as a function of primary energy. Results from our simulation 
for air spheres of radius R = 1 cm and R = 10 cm are shown (solid lines) 
for comparison with the energy deposition values reported by other authors 
(symbols and broken lines). See text for details. 
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for 28.5 GeV electrons crossing a 1.67 cm air gap using EGS4 which is closer 
to ours. In next section a detailed comparison is shown. 

The MACFLY Collaboration [13], using GEANT4, reported a (dE/dx) dep 
curve obtained for a cylindrical chamber, 15 cm in diameter and 28.8 cm 
long. Its results above 2 MeV are in fair agreement with our predictions 
for a sphere of radius 10 cm at atmospheric pressure, except for the highest 
energies (> 100 MeV), where our values are up to 6% larger. At lower 
energies the values are oddly much larger than expected (see also the next 
section for more details). 

Finally, the AIRFLY Collaboration, using GEANT4, has provided us 
with the results of a dedicated simulation for our simple geometry (sphere of 
radius 1 cm and primary electrons are forced to interact in its center) [30J, 
which are fully compatible with ours within statistical fluctuations. 

3 Detailed simulation of fluorescence yield 
experiments 

In this section, a simulation including the particular geometry for several 
experiments as well as the propagation of beam electrons will be presented. 
That allowed the calculation of the energy deposited in the corresponding 
observation volume. When possible, the average geometrical acceptance for 
fluorescence detection has been evaluated and compared with that assumed 
by the authors of the experiment. For that purpose, the solid angle sub- 
tended by the optical system's entrance pupil is averaged along the electron 
trajectories, weighting with the amount of light emitted at each point. Note 
that the effect on the detection efficiency due to the spread of electrons in 
the interaction volume is presumably small, as it results from the average 
of a distribution with axial symmetry around the beam axis. That is, the 
enhanced fluorescence contribution from electrons approaching the fluores- 
cence detector is nearly compensated by the reduced one from those moving 
in the opposite direction. In addition, in those cases where the electron- 
photon coincidence technique is used (e.g., with a radioactive source), the 
detection of beam electrons is required in the simulation, and thus, primaries 
not firing the detector are rejected. This may affect the energy spectrum of 
contributing electrons, and so, the determination of the deposited energy. 
We have performed this detailed simulation for the experiments of Nagano 
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et al. HaHSl, AirLight [IB], FLASH-2008 [20] and MACFLY [13]. Experi- 
ments of Kakimoto et al. [H] and Lefeuvre et al. [T7j will be also discussed 
although a dedicated simulation was not carried out. 

3.1 Nagano's experiment 

For Nagano's setup [151 HB], we made a detailed simulation in order to de- 
termine the effect of several approximations applied in this experiment and 
then evaluate a correction to their results. The collisional chamber of this 
experiment is shown in figure [2j where the geometry implemented in this 
work is emphasized. They used the electron-photon coincidence technique, 
so this feature has been also implemented in our simulation. Electrons fire a 
scintillator at a distance of about 16 cm from the 90 Sr source. However, only 
the 5 cm gap between the collimator and the scintillator enclosure is in the 
field of view of the photon detectors. In this setup the observation volume 
corresponds to the chamber volume excluding the collimation channel and 
the scintillator enclosure. 

As input to the simulation, we have used a truncated gaussian fit (solid 
line in figure [3]) to the measured 90 Sr spectrum of [15], with mean energy of 
0.85 MeV, end point at about 2 MeV and a threshold of 0.3 MeV due to the 
material covering the radioactive source. According to our simulation, most 
beam electrons are scattered away at atmospheric pressure and so are not 
detected. This occurs preferably for low-energy primaries, and therefore the 
spectrum of contributing electrons is shifted to higher energies, with a mean 
value of 1.11 MeV at atmospheric pressure (see figure [3]). 

The approximations made by Nagano et al. are the following. Firstly, 
they assumed that all the fluorescence light is emitted from the beam it- 
self. However, as explained above, a non-negligible fraction of fluorescence is 
produced by secondary electrons at some distance, and thus, the average ge- 
ometrical acceptance Q may be different from that which would be obtained 
if all the light were generated in the beam f2 beam . Secondly, they calculated 
the e value assuming the photons are emitted in a length equal to the gap 
distance Ax gap = 5 cm. However, due to the fluctuating trajectories of beam 
electrons, the mean path length Ax should be somewhat larger. Finally, 
for the calculation of the fluorescence yield from ([2]), they assumed that the 
energy deposited in the observation volume equals the collisional energy loss 
(d-E/dx)i oss at 0.85 MeV as given by the Bethe-Bloch formula. As discussed 
above, the deposited energy is expected to be smaller than the total energy 
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photomultiplier tube 




Figure 2: Schematic drawing of the chamber of Nagano's experiment (taken 
from [IS])- Electrons from a 90 Sr source are beamed and detected by a 
scintillator counter. Between the collimator and the scintillator enclosure 
there is a 5 cm gap in the field of view of the photon detectors. The thick 
contour represents the geometry implemented in our simulation. The energy 
deposition obtained from our simulation including these geometrical details 
is fully compatible with the prediction for a sphere of 5 cm radius (circle). 
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Figure 3: Energy spectrum used as input to the simulation of Nagano's 
experiment (solid line), with a mean energy of 0.85 MeV. According to our 
simulation, the energy distribution of electrons reaching the scintillator (sym- 
bols) is biased to larger energies, resulting in a mean value of about 1.11 MeV 
at atmospheric pressure. Both spectra are normalized to the same area. 
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loss. 

Therefore, the fluorescence yield value YNag reported by Nagano et ai, 
should be corrected by three factors accounting for the above described ap- 
proximations 

^bcam Ax ga p (dE/dx) loss 

r " rNas n Ax (dE/dx) Aep ' 1 ' 

From our detailed simulation, we have obtained the average geometrical 
acceptance relative to that assumed in [151 IB], that is, the f2/f2beam ratio in 
the above expression. For this calculation, we have taken into account that 
some regions in the observation volume are hidden to the optical system. We 
have found that this correction to the efficiency is very small, as expected, 
resulting Q to be about 1% lower than f2 beam at atmospheric pressure. The 
correction is even smaller at lower pressure. On the other hand, according 
to our simulation, the mean path length of contributing electrons at atmo- 
spheric pressure is about 1% longer than the gap length. Therefore, these 
small corrections happen to compensate each other in equation ^ for this 
experiment. In fact, the Nagano's results on e (photons/meter) turn out to 
be unaffected by these approximations. 

The main correction to the fluorescence yield in equation ^ arises from 
the approximation in the energy deposition. From our simulation, we have 
obtained the average energy deposited per electron and unit path length 
(d-E/dx)dep as the ratio between the integrated energy deposition and the 
total path length of beam electrons within the observation volume. Simula- 
tion data were accumulated up to reach a statistical uncertainty of less than 
1% in the {dE/dx)d ep value. As a result of this correction, the fluorescence 
yield at atmospheric pressure from [151 US] should be increased by 7%. 

The result of the deposited energy from our simulation of the Nagano's 
experiment is fully compatible with the predictions from the generic simu- 
lation for a simple geometry, assuming a sphere of radius R = 5 cm (see 
figure [2]) and an electron energy of around 1 MeV. As already mentioned, 
deposited energy is very weakly dependent on R. In addition, at these en- 
ergies around the minimum of the Bethe-Bloch formula, the result is nearly 
independent of E. Therefore, for experiments similar to that of Nagano, the 
generic simulation may provide reliable results on the energy deposited in 
the observation volume. 
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3.2 The AirLight experiment 

The experimental technique used in AirLight [18] is similar to the one of 
Nagano, but with a collisional chamber larger by about a factor of two. 
Electrons from a 90 Sr source are collimated to form a narrow beam of 5 mm 
diameter. Between the collimator and the electron detector (a scintillator 
of 4 cm diameter) there is a 10 cm gap within the field of view of seven 
photomultipliers. In [18] a detailed simulation of the experiment is carried 
out using GEANT4, including electron scattering by elements of the cham- 
ber (e.g., collimator walls, scintillator). Both beam electrons and secondary 
ones are thoroughly tracked allowing the authors the determination of the 
geometrical acceptance as well as the energy deposited in the observation 
volume. 

For this experiment we have carried out a simulation including primary 
propagation and the geometrical details of the setup. Electron scattering by 
other elements of the chamber apart from air is not included. The aim of 
our simulation has been the calculation of the integrated energy deposition 
as a function of the electron energy to be compared with the results found by 
AirLight. While at low pressure we have obtained results in fair agreement 
with those reported in [18] (within 5% at 50 hPa), at atmospheric pressure 
our simulation predicts a larger energy deposition. For instance at 800 hPa, 
deviations range from about 5% at low electron energy up to 20% at 2000 keV 
(see figure [4]). The origin of this discrepancy is still unclear. It should be 
indicated that electron scattering by components of the collision chamber 
(included in their simulation, but not in ours) would result indeed in an 
increase of the energy deposition, as pointed out in [TH] . 

From our simulation we have also calculated the integrated value of the 
collisional energy loss. At atmospheric pressure we obtain results somewhat 
larger than those reported by [TH], in particular at low electron energy (see 
figure |4|. The energy loss computed by AirLight seems to have been calcu- 
lated as the product of the Bethe-Bloch stopping power (dE / dx)\ oss [2H] times 
the 10 cm gap length, and thus, assuming that the stopping power is con- 
stant along the electron path. However, as already mentioned for Nagano's 
experiment, due to the fluctuating trajectories of primaries, their mean path 
length is expected to be somewhat larger than the gap length. Also, at low 
electron energy, the average stopping power should be larger than the initial 
value, because primaries lose a significant fraction of their energy in their 
path. 
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Figure 4: Integrated energy deposition in the chamber of the AirLight exper- 
iment as a function of the electron energy detected in the scintillator in the 
200 - 2000 keV range p]. Results at 800 hPa obtained by AirLight using 
GEANT4 (squares) are compared with ours (circles connected by lines). In- 
tegrated energy loss from our simulation (dotted line) is compared with that 
calculated by AirLight (solid line). See text for more details. 
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Obviously, the energy deposition must be smaller than the energy loss, 
and their difference is expected to diminish at low E values. This behavior 
is found in our simulation as well as in that of [18]. However, at 800 hPa, 
the energy loss from [18] is smaller than our prediction of the energy de- 
position for electron energies below 400 keV (figure ^ showing an evident 
inconsistency between both calculations. 

Assuming a mean electron energy of 1 MeV, our detailed simulation pre- 
dicts an average deposited energy around 10% higher than that of [18] at 
800 hPa, and thus, if our results were correct, their fluorescence yield value 
at this pressure should be reduced in principle by the same factor. However, 
AirLight reports Y® v , results (together with a set of quenching parameters) 
by averaging measurements carried out at several pressure^] In fact the 
average Y® v , values are nearly equal to those obtained at 400 hPa. At this 
pressure and an electron energy of 1 MeV, the difference between our energy 
deposition and that of AirLight is about 7%. Therefore, according to our 
simulation, results on fluorescence yield from [T8] should be reduced by an 
overall correction factor of 7% irrespective of pressure. 

3.3 The FLASH experiment 

In this setup [20], an intense beam of 28.5 GeV electrons from a linear accel- 
erator crosses a cylindrical collision chamber, 25 cm in length and 15 cm in 
diameter. The beam is coaxial to a pair of 1.6 cm diameter cylinders with 
a 1.67 cm gap between them, which defines the beam path length as seen 
by two PMTs. Each detector is mounted at the end of a long channel per- 
pendicular to the beam axis and provided with baffles to suppress unwanted 
scattered light. This channel terminates, at 45 cm from the beam, in a 1.2 cm 
diameter window that forms the entrance pupil of the optical system. In the 
opposite side, there is a shorter channel holding a LED used to monitor the 
PMT gain stability, being the extreme of this channel at about 27 cm from 
the beam (72 cm from the window of the PMT). The absolute calibration of 
the optical system was carried out by comparison with the measurement of 
Rayleigh-scattered light from a narrow beam of a nitrogen laser. As pointed 
out above, the FLASH Collaboration used an EGS4 simulation to obtain 
the mean energy deposited per electron in the observation volume of their 
experiment [20J. Their simulation also provided the average geometrical ac- 

6 AirLight found a slight, although unexpected, pressure dependence of Y° v ,. 
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ceptance for fluorescence detection relative to that for the calibration laser 
beam. 

We have carried out a detailed simulation of this experiment in order 
to evaluate the energy deposited in the chamber as well as the average ge- 
ometrical acceptance for comparison with their results. In our algorithm, 
the observation volume has been approximated by a rectangular prism of 
dimensions 1.60 cm x 1.67 cmx72 cm. 

The deposited energy is not explicitly given in [20], but it can be inferred 
from the ratio between the reported results on e and Y at pressures ranging 
from 67 to 1013 hPa. In figure |5j the resulting values are shown as a func- 
tion of pressure together with those from our detailed simulation^, where 
statistical uncertainties were reduced to less than 0.3%. At high pressures 
(> 500 hPa), both results are in good agreement (within 2%). However, the 
pressure dependency is different, resulting in a deviation of 6% at 67 hPa. 
This discrepancy could be due to a different treatment of the density correc- 
tion, which at such a high electron energy is very important, lowering the 
inelastic cross sections by about 40% at atmospheric pressure. As a simple 
test, we have recalculated the deposited energy using in our simulation the 
density correction for 1013 hPa in the whole pressure range. The resulting 
energy deposition, also presented in figure [5| shows a behavior similar to that 

of nu. 

In regard with the geometrical acceptance, FLASH [20J found that the 
efficiency of their optical system for fluorescence detection is (3.2 ± 0.25)% 
lower than that for the calibration laser beam, due to the wide-spread energy 
deposition from the electron beam. However, our simulation predicts that 
this effect is negligible in the whole pressure range for this setup. 

Notice that according to our simulation both corrections, i.e., deposited 
energy and optical efficiency, nearly cancel and therefore the fluorescence 
yield at 1013 hPa from [20] would be almost unaffected (see table [l] in sec- 
tion |5|). 

3.4 The MACFLY experiment 

In this experiment [13], the collision chamber is a cylinder 28.8 cm long and 
15 cm in diameter with a mirrored inner surface in order to increase the ge- 

7 The results of our detailed simulation are in excellent agreement (within 1%) with 
that obtained from the generic simulation for a sphere of 1.67 cm radius. 
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Figure 5: Energy deposited per electron and unit mass thickness as a func- 
tion of pressure for the FLASH experiment [20] ■ Results from [20] (squares 
connected by lines) and from our detailed simulation (solid line) are in good 
agreement at atmospheric pressure, although the pressure dependency is dif- 
ferent. For comparison, results from our simulation applying the density 
correction at 1013 hPa in the whole pressure range are also shown (broken 
line). 
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ometrical acceptance. They used the coincidence technique, with individual 
electrons either from a linear accelerator (20 and 50 GeV) or from a 90 Sr 
source of mean energy 1.5 MeV. As mentioned above, MACFLY calculated 
the energy deposited per electron and per unit mass thickness as a function 
of E using GEANT4 (blue dotted line in figure Q. 

We have carried out a detailed simulation to obtain the energy deposition. 
Primaries colliding the walls of the chamber or being stopped inside, and so 
not reaching the trigger system, are rejected in the simulation. The average 
energy deposited per electron and unit mass thickness has been obtained from 
the ratio between the integrated energy deposition and the total path length 
of triggering primaries. Statistical uncertainties are reduced below 0.5%. The 
deposited energy obtained from our simulation is slightly lower (2%) than the 
value reported by MACFLY at 1.5 MeV and larger by about 6% at 20 and 
50 GeV, and therefore, according to our calculations the fluorescence yield 
of MACFLY should be increased by 2% at 1.5 MeV and decreased by 6% at 
20 and 50 GeV. 

The above described detailed simulation is in good agreement (better than 
2%) with the generic simulation for R = 10 cm at energies above 0.5 MeV 
(see figure [If. Below 0.5 MeV, primary electrons lose a significant fraction of 
their energy in the chamber (e.g., electrons with energy lower than 0.2 MeV 
are fully stopped), and consequently the (d£'/da:)dep values resulting from our 
detailed simulation are somewhat larger (~ 5%) than those obtained from the 
generic one. However this effect cannot justify the very large discrepancies 
with the MACFLY results at low electron energy shown in figure [Tj 

3.5 Other experiments 

Apart from those discussed in the above sections other well-know experi- 
ments have been carried out, in particular those of Kakimoto et al. [H] and 
Lefeuvre et al. [ITj . For these we have not performed a detailed simulation 
but we will use the result of our generic simulation. 

In the Kakimoto's experiment electrons with a mean energy of 1.4 MeV 
from a 90 Sr source as well as those from an electron synchrotron with energies 
of 300, 650 and 1000 MeV were used. Fluorescence light was produced and 
observed in a 10 cm gap. For the determination of the fluorescence yield, 
they assumed full energy deposition in the beam axis. According to our 
generic simulation, the deposited energy inside a sphere of 10 cm radius 
is about 6%, 25%, 28% and 29% lower than the energy loss for 1.4, 300, 
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650 and 1000 MeV electrons, respectively. Therefore, the fluorescence yield 
values at the above energies should be increased correspondingly by these 
factors. Results obtained from applying these corrections to measurements of 
Kakimoto et al. are consistent with those obtained in [31] using an alternative 
procedure. 

In the experiment of Lefeuvre et al. |17j . electrons with mean energies 
of 1.1 MeV (or 1.5 MeV when applied a higher energy threshold to the 
electron detector) from a very intense 90 Sr source (370 MBq) were used. In 
order to shield the fluorescence detector from the large amount of X-radiation 
produced at such a high electron rate, both the radioactive source and the 
scintillator counter were placed inside a solid lead cylinder, 10 cm in diame- 
ter. Electrons trajectories were confined within a cone dug in the lead and 
coaxial to the cylinder, being the 90 Sr source at the apex and the scintillator 
on the base. A 4 cm diameter hole was bored through the lead cylinder, 
allowing fluorescence to be measured from both sides. In order to obtain the 
average geometrical acceptance and the mean path length of electrons, they 
performed a GEANT Monte Carlo simulation including electron scattering 
by the lead walls. They found that this feature has an important role in 
this experiment, since half of primaries reach the scintillator after hitting 
the cone walls. From their simulation, they also found that a very small 
fraction of secondary electrons have enough energy to escape the observation 
region at atmospheric pressure. In particular, they obtained that only 1% 
of secondaries have energies above 5 keV, and 33% of them are lost in the 
lead before having produced any fluorescence. From these simulation results, 
they inferred that the deposited energy nearly equals the total energy loss 
and a minor correction of 0.4% was applied. However, as discussed above, a 
significant fraction of energy is deposited by these few energetic secondaries, 
and thus, the correction is expected to be larger. 

From our generic simulation assuming a sphere of 4 cm radius for Lefeu- 
vre's experiment, we have obtained that the deposited energy is 9% and 
10% lower than the total energy loss at mean electron energies of 1.1 and 
1.5 MeV, respectively. However, electron scattering by the lead shield of [T7| 
is not included in our simulation, which assumes that secondaries reaching 
the walls are absorbed, and thus, it may underestimate the deposited energy 
in this case. We have evaluated this effect in order to estimate a correction 
to our simulation results for this experiment. Using CASIN02.42 [32J, we 
have found that, for a 60° incident angle and energies in the 1 keV - 1 MeV 
range, about 60% of electrons are backscattered by lead, having lost almost 
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half of their energy on the average. Assuming that backscattered electrons 
deposit all their remainder energy in the air volume, about 30% of the energy 
of electrons reaching the lead walls is recovered, and thus, should be added to 
the total energy deposition. As a consequence, the fluorescence yield values 
at 1.1 and 1.5 MeV from [T7] should be increased by about 7% according to 
our calculations. 



4 Wavelength interval 

The wavelength normalization can be carried out if the relative intensities of 
the molecular bands within the experimental interval are known. As is well 
known, relative intensities of bands sharing the same upper level are given 
by the relative Einstein coefficients B vv i = A vv r / ^2 V , A vv i . As described 
below fluorescence intensities of bands with different upper levels can be also 
calculated from molecular parameters. 

In the absence of other population mechanism^] apart from direct exci- 
tation of the vibrational level v from the ground state X, the band inten- 
sities of a given system are proportional to the excitation cross sections of 
the corresponding upper levels. At high electron energy (i.e., in the Born- 
Bethe high-energy limit), these cross sections are known to be approximately 
proportional to the Franck-Condon factors qx^v, which correspond to the 
probability ratios for optical transitions. It has been found [10] that this 
approximation can be also applied in our case, in spite of the fact that most 
fluorescence light is generated by low-energy secondaries. Under this approx- 
imation, for a given band system, the relative intensities I vv > with respect to 
a reference one Iqq (i.e., the most intense band of the system, located at 
337 nm for the 2P system and at 391 nm for the IN system) are given by 

Qx->vB vv > 1 + P/Pq , _v 

loo ~ gx^o^oo 1 + P/P v ' 1 J 

At atmospheric conditions, where P » P,| can be assumed, the above 
expression becomes 



Iw' Qx— >vB vv i P v 



ho Qx^qBooPq 



(8) 



8 



As indicated in section 



2.1 



the effect of vibrational relaxation of the 2P upper levels 



can be considered as included in the P' v quenching parameters. 
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where any pressure or temperature dependence is cancelled, since the same 
approximate y/T growing behavior can be assumed for all P' v parameters. 

Relationship (J8|), and hence the applicability of the Franck-Condon fac- 
tors, has been verified using experimental data on relative intensities. In fig- 
ure 2 of ref. [10J we plotted the relative intensities measured by AIRFLY [25] 
against q x ^ v B vv iP' v for the 2P system, using the tabulated values of qx^v 
and B vv i from [33], EI] and the experimental characteristic pressures for dry 
air reported in [25J. A remarkable linearity was found, in particular for the 
bands with upper vibrational level v — 0, 1, contributing to almost 90% of 
the total 2P fluorescence. A similar analysis was previously performed in [H] 
using the relative intensities and P' v values measured by Nagano et al. [16J. 
Although linearity was also observed, significant fluctuations were found very 
likely due to much larger experimental uncertainties than those of |25j . 

Expression (|8]) can be extended for molecular bands belonging to differ- 
ent systems including the ratio between the respective effective cross sec- 
tions cXgy St , defined as a^, = al^ st qx^ v B vv i [ID]- Unfortunately, whereas the 
Franck-Condon factors and the Einstein coefficients are known very precisely, 
the available data on <rfy St and P' v are less reliable. Effective cross sections de- 
pend on molecular parameters with significant uncertainties and are geometry 
dependent, and the available experimental data on characteristic pressures 
show large discrepancies. However, as already pointed out in [TU] , the relative 
intensities measured by different authors, which are proportional to a^ st P' v) 
are in good agreement. For instance, for the /337//391 ratio at atmospheric 
pressure, values of 3.0 ± 0.7, 3.7 ± 0.4, 3.4 ± 0.3 and 3.57 ± 0.13 result from 
data from [35], [15], [16] and [25J, respectively. Even a very consistent 3.5 
value was measured by Bunner pQ for 4 MeV a-particles as well. Therefore, 
the value of the <7^ t P^ products from experimental relative intensities is very 
reliable. 

In order to compare absolute results on fluorescence yield from different 
authors, we will reduce values measured in a given wavelength interval A A 
to the yield at the 337 nm band. Experimental results are then normalized 
multiplying by the ratio between the relative intensity of the 337 nm band 
and the sum of all the bands within A A. This /337//AA ratio can be thus 
written in terms of the empirically reliable /337//391 ratio as 

-^337 
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where relative intensities at atmospheric conditions follow ([8]) separately for 
each band system. Note that, due to the 2P dominance, even an error of a 
factor 2 in the I337//391 ratio results in an uncertainty below 10% in /337//AA 
for the spectral ranges of the experiments considered here. 

For this wavelength reduction, we use accurate experimental relative in- 
tensities from [22] in the 284 - 429 nm spectral range, including some bands 
of the Gaydon- Herman (GH) system [25] (1.6% total contribution). For 
instance, for the whole wavelength interval of [25] a ratio 1337/1(284-429] = 
0.2569 ± 0.0022 is obtained, where the uncertainties on the relative intensi- 
ties are treated as uncorrelated as the authors suggest. Note that this value 
is slightly lower than those listed in table 1 of ref. JTU], where intensities of 
bands of the GH system were neglected. In addition, equations ^ and ^ 
allow the determination of /337//AA for a wavelength interval including some 
of the 2P and IN weak bands outside the spectral range of |25j . 

5 Results of the fluorescence yield 
comparison 

In table [T] absolute fluorescence yield values at atmospheric pressure reported 
by several authors (first column) are put in comparison. The wavelength in- 
terval, pressure, temperature and electron energy at which absolute values 
are reported are shown in columns 2-5. The experimental results, in units 
of either photons/MeV or photons/m, and the corresponding uncertainties 
provided by the authors are quoted in columns 6 and 7. In order to get 
values normalized to that of the 337 nm band we have applied the procedure 
shown in section |4j The corresponding normalization factors are displayed 
in column 8. The last column shows the fluorescence yield for the 337 nm 
band in units of photons/ Me V. For an easy comparison, all results are given 
for 800 hPa and 293 K (dry air). For this normalization we have used ex- 
pression (3j assuming P 3> P 337 and P 337 oc s/T. The numbers in the left 
side of this column correspond to the Y337 value obtained using the energy 
deposition and geometrical factors assumed (or calculated) by the authors. 
Finally, in the right side (in bold) we show the results we have obtained 
applying the various corrections discussed in section |3j 



As explained in section 2.3, our simulation gives us both deposited energy 



and fluorescence emission and thus we can provide a theoretical result of 
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the fluorescence yield [36J. The updated calculations of this work lead us 
to a Y" 3 3 7 value of 390 ph/MeV, nearly independent of the electron energy 
above 0.1 MeV. This theoretical result in combination with the characteristic 
pressure for the 337 nm band in dry air [25] gives rise to a Y337 value of 
7.5 ph/MeV at 800 hPa and 293 K, which is close to the experimental ones 
shown in the table. However, as pointed out above, uncertainties in the 
simulated fluorescence emission are large, and thus, this theoretical result 
cannot be considered reference one. 

6 Conclusions 

A direct comparison of the absolute measurements of the air fluorescence 
yield available in the literature is not possible since very often the reported 
data correspond to different wavelength intervals and the magnitude is ex- 
pressed in different units. The conversion factors for such comparison depend 
on the particular features of the experimental setup. In this work we have de- 
scribed a procedure for the normalization of the various experimental results 
to a common unit and wavelength interval. 

In the first place, fluorescence yields have been reduced to that corre- 
sponding to the most intense band at 337 nm using accurate experimental 
data on intensity ratios in full agreement with theoretical relationships. Ab- 
solute values can be expressed in photons/Me V units as far as the energy 
deposited in the field of view of the setup is known. In this work we have 
used an upgraded simulation algorithm to calculate the energy deposition of 
electrons in a wide energy range. A generic algorithm has provided results 
as a function of the pressure and the overall size of the collision chamber. In 
addition, for several experiments, we have performed a detailed simulation 
(in full agreement with the predictions of the generic one) including in some 
cases the calculation of geometrical factors. The results of these simulations 
have been compared with those reported by the authors of the corresponding 
experiments and eventual disagreements have been discussed. 

As a result of this work a table of the absolute fluorescence yields normal- 
ized to 800 hPa and 293 K have been presented. According to our simulation, 
a non-negligible correction has to be applied to those results obtained un- 
der the assumption that the energy lost by the primary electrons is fully 
deposited in the field of view of the optical system (i.e., Kakimoto et al. [14], 
Nagano et al. [T5JEES]) an d Lefeuvre et al. [H]). In particular the Nagano's 
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fluorescence yield should be increased by about 7%. Possible implications 
of these corrections to the Y values in the determination of the end of the 
energy spectrum of cosmic rays have been discussed recently in [37J. 

On the other hand, our calculations of energy deposition are in general 
in reasonable agreement with those from experiments that carried out a de- 
tailed simulation of their setup. Even though, in some cases we have found 
non-negligible discrepancies and correction factors have been proposed. Ex- 
periments claiming high accuracy should be very careful in the evaluation of 
the error contribution from the deposited energy. 

From the above table, and assuming the corrections proposed in this 
work, most fluorescence yield measurements at those reference conditions 
are contained within the 6-7 ph/MeV interval. Notice, however that the 
fluorescence yield reported by MACFLY, as well the preliminary result of 
AIRFLY, are smaller by more than 10%. 
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